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Abstract 

We extend the framework of Nambu-Hamiltonian Mechanics to include dissipation in i? 3 
phase space. We demonstrate that it accommodates the phase space dynamics of low dimen- 
sional dissipative systems such as the much studied Lorenz and Rossler Strange attractors, as 
well as the more recent constructions of Chen and Leipnik-Newton. The rotational, volume 
preserving part of the flow preserves in time a family of two intersecting surfaces, the so 
called Nambu Hamiltonians. They foliate the entire phase space and are, in turn, deformed 
in time by Dissipation which represents their irrotational part of the flow. It is given by the 
gradient of a scalar function and is responsible for the emergence of the Strange Attractors. 

Based on our recent work on Quantum Nambu Mechanics, we provide an explicit quan- 
tization of the Lorenz attractor through the introduction of Non-commutative phase space 
coordinates as Hermitian N x N matrices in R 3 . They satisfy the commutation relations 
induced by one of the two Nambu Hamiltonians, the second one generating a unique time 
evolution. Dissipation is incorporated quantum mechanically in a self-consistent way having 
the correct classical limit without the introduction of external degrees of freedom. Due to its 
volume phase space contraction it violates the quantum commutation relations. We demon- 
strate that the Heisenberg-Nambu evolution equations for the Quantum Lorenz system give 
rise to an attracting ellipsoid in the 3N 2 dimensional phase space. 
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1 Introduction 

Dissipative dynamical systems, with a low dimensional phase space, present an important 
class of simple non-linear physical systems with intrinsic complex behavior (homoclinic bi- 
furcations, period doubling, onset of chaos, turbulence,...), which generated intense experi- 
mental, theoretical and numerical work in the last few decades [TJ [2j [3J S] - 

In this work we shall be interested to study dissipative dynamical systems with a 3- 
dimensional phase space from the perspective of Nambu-Hamiltonian Mechanics(NHM) [5, 6j. 
It represents a generalization of Classical Hamiltonian Mechanics, which is mostly appropri- 
ate for the study of phase-space volume preserving flows ( Liouville's theorem). To study 
dissipative systems in R 3 we generalize NHM by splitting the flow vector field into its rota- 
tional(solenoidal) and irrotational components, which we identify them to be its integrable, 
non-dissipative as well as dissipative parts respectively. 

We apply this idea to the most famous examples of Lorenz [7] and Rossler [8] chaotic 
attractors which represent the prototype models for the onset of turbulence. EU EJ H2j . 
We present also an Euler Top deformation of the Lorenz system |13| [T4] which will be a 
reference system for the quantization of the chaotic attractors. Following our recent work on 
the Quantization of NHM |15| IT6] we provide a Matrix model for N-interacting Lorenz ( or 
Rossler) attractors. It replaces the classical three phase space coordinates by three N x N 
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Hermitian matrices and the evolution equations with a symmetrically ordered matrix system 
(Weyl Quantization) [17 . As a consequence we obtain a chaotic dynamical system of 3iV 2 
real dimensional phase space for any integer N. Finally we demonstrate that, similarly to the 
Lorenz system, there is an attracting ellipsoid in the 3iV 2 -dim phase space, which provides a 
compact hypersurface for the localization of the higher dimensional attractors. 
The plan of our paper is as follows: 

In sect. 2 we introduce the framework of Dissipative Nambu-Hamiltonian Mechanics in 
3-dim phase space i? 3 . Its Dynamics decomposes into a Non-dissipative integrable sector, 
which is parametrized in terms of two intersecting surfaces, the Nambu Hamiltonians (or 
Clebsch-Monge potentials of hydrodynamics) defined at each point of the phase space flow. 
One of them defines (together with the initial conditions) a two dimensional phase space 
embedded in i? 3 , while the other defines the time evolution of the system, on this curved 
2d-phase space. Their intersection defines the trajectory of the system and the flow vector 
field. The dissipative part of the flow is the gradient of a scalar function. 

In sect. 3 we apply the general framework to the Lorenz system. We isolate a conservative 
sector parametrized by two intersecting surfaces: a circular and a parabolic cylinder . They 
account, amazingly, for the double scroll topology of the full " butterfly" attractor. As 
a consequence of the linear character of dissipation the trajectories of the full system get 
organized around the trajectories of the exactly integrable part. 

The Nambu-Hamiltonian form of the integrable part, given also in terms of Poisson brack- 
ets in ii 3 , defines, on one of the surfaces, a two dimensional phase space (the parabolic cylin- 
der) with a standard Poisson structure, of a quartic anharmonic oscillator of single or double 
potential well, depending on the initial conditions. 

We present also an Euler top deformation of the Lorenz system . It is generated by a 
conservative sector given by the intersection of a cylinder with a paraboloid. 

In sect. 4 we apply our method method of analysis to the well known case of the Rossler 
attractor. In this case the non-linear dissipative part distorts completely the orbits of the 
non-dissipative sector given by the intersection of a cylinder with an helicoid. 

In sect. 5 we describe the basic steps of our Quantization prescription for Dissipative 
Nambu-Hamiltonian Mechanics. 

In sect. 6 we construct in detail the Matrix Model for N interacting attractors as an illus- 
tration of our Quantization method. Moreover we demonstrate the existence of an attracting 
ellipsoid in 3iV 2 dim. phase space. 

Finally in sect. 7 we conclude our work with the basic results and a discussion of 
the relevance of our proposal to the formal quantization problem of dissipation of chaotic 
systems [Eldn]. 
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2 Dissipative Systems in R 3 and the Dynamics of 
Intersecting Surfaces 



Nambu-Hamiltonian mechanics is a specific generalization of classical Hamiltonian mechan- 
ics, where the invariance group of canonical symplectic transformations of the Hamiltonian 
evolution equations in 2n dimensional phase space is extended to the more general volume 
preserving transformation group SDiff{M) with an arbitrary phase space manifold M of 
any dimension d = dim(.M ) . 

In the present paper we will work with the case of A4 = R 3 , a three dimensional flat 
phase space manifold. Nevertheless our results could be generalized to curved manifolds of 
any dimension |6j IT5] . 

Nambu-Hamiltonian mechanics of a particular dynamical system in R 3 is defined once 
two scalar functions H{ £ C°°{R 3 ),i = 1, 2, the generalized Hamiltonians (5j |6] are provided. 
The evolution equations are: 

x { = {x i ,H 1 ,H 2 } i = 1,2,3 (2.1) 

where the Nambu 3-bracket, a generalization of Poisson bracket in Hamiltonian mechanics, 
is defined as 

{f,g,h} = ^ k d l f&gd k h i,j,k= 1,2,3 V/, g, h e C°°(R 3 ) (2.2) 
Any local coordinate transformation 

x i - y l = y\x) z = 1,2,3 (2.3) 

which preserves the volume of phase space 

d6t (IS") =1 ^ x = (x\x 2 ,x 3 ) e R 3 (2.4) 

leaves invariant the 3-bracket and therefor it is a symmetry of Nambu-Mechanics. Except 
for the linearity and antisymmetry of the bracket with respect to all of its arguments it also 
satisfies an important identity, the so called "Fundamental identity" [FI] [Bj. Let us define the 
operator 

C( hlM )f = {h!,h 2 J}, V hx,h 2 ,f GC°°(R 3 ) (2.5) 
Then it can be shown that C satisfies the commutation relation 

[£(fci,h 2 ) > £(/i3,/u)] = £{hi,h 2 ,h 3 },h4 + £-hz,{hiMM} V /ij £ C°°(i? 3 ), i = 1, 2, 3, 4 (2.6) 

In |15| it was shown that rel(2.6) are the commutation relations of the Lie algebras of Volume 
preserving diffeomorphism Group SDiff(i? 3 ). The argument goes as follows: First using the 
operator J0m lt jj 2 \ the evolution equations , rel (2.1), can be written as 

x* = C Hl ,H 2 x l i = 1,2,3 (2.7) 
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Given any function / 6 C°°(-R 3 ) of the phase space coordinates we obtain the Liouville eq. 

/ = C Hl , H2 f = {/,#!, #2 } (2.8) 

This equation demonstrates that Hi,H 2 are constants of the motion and the orbit is given 
by the intersection of the surfaces H° = H°, H 2 = H 2 where 

H° x = H^xl), H° = H 2 (xl) (2.9) 

with x l Q specifying the initial conditions at t = 0, x l Q = x l (0). The formal integration of 
rel.(2.9) is: 

/(s*(t)) = e tC ^ H °^ f(xl) (2.10) 
The evolution eq.(2.1) has a flow vector field: 

v *( x ) = e ljk d j H 1 d k H 2 i,j,k = 1,2,3 (2.11) 

which is volume preserving 

dV = (2.12) 

The reverse is also true. Any flow vector field of zero divergence v l G C°°(R i ), i = 1,2,3 
can be represented at least locally by two Clebsch-Monge potentials H\ , H 2 |15| 120"! |2~H [22] 

y i = e ijk d j Hld k H2 (2.13) 

Prom the representation rel.(2.13) and the commutation relation for the Lie algebra of 
SDiff(i? 3 ) : 

[£ u ,£v } = C w (2.14) 
where C u = u % d % (similarly for (v, w) with d % u l = d l v t = and 

w i = u-W-v-du* (2.15) 

follows that d l w l = and hence the commutation relations (2.6) follow. We will name 
the phase-space volume preserving flows " Non-dissipative" while the non-conserving ones 
" Dissipative" . 

The general flow vector fields will have regions D + C B? with 

<9V(x) > x£D + (2.16) 
where any volume element is expanded by the flow and regions L>_ C R 3 , D + n D- = with 

3V(x) < VieD- (2.17) 
where the volume elements are contracted. 



The standard local stability analysis of the flow, involves the determination of the critical 
points x* (equilibrium points) where v l {x ) = 0, i = 1,2,3, as well as the examination of 
the eigenvalues of the fluctuation matrix (6>V),£,j = 1,2,3) at the critical points which 
characterize the phase space portrait of the flow. 

Under variations of external parameters (environment) the location and stability character 
of the critical points vary. There arises possibilities of Hopf bifurcations, transition to chaos, 
turbulence etc [H E3 G9 15 G3 . 

In order to bring into play the Nambu-Hamiltonian Mechanics framework we firstly gen- 
eralize it for general dissipative systems in B? phase space. We will decompose any flow 
vector field into its rotational and irrotational parts. 

v i = e ijk d j A k + d i D i,j,k = l,2,3 (2.18) 

with the vector potential determined up to a local gauge transformation: 

A* -»■ A i + a* A i = 1,2,3 (2.19) 

and 

0V = V 2 D (2.20) 

Given d l v l and appropriate BC we can , in principle, determine the dissipation function D. 

In general D is determined up to harmonic functions on R s . Since the rotational com- 
ponent is diver genceless, we shall call it the Non- dissipative part of the flow (ND), vnd 
while the irrotational component we will identify it with the Dissipative part(D), vp. Once 
the vector potential A 1 has been determined we can choose the Clebsch-Monge gauge where 
[CM]: 

A i = H x d i H 2 (2.21) 

utilizing two "generalized Hamiltonians" or equivalently Clebsch-Monge flow potentials. The 
flow equations take the form 

x i = {x i ,H 1 ,H 2 } + d { D t = 1,2,3 (2.22) 

and in vector form 

'x = ViTi x VH 2 + VD (2.23) 
We notice that given a pair of functions Hi , H 2 such that 

v ND = V#i x VH 2 (2.24) 

any transformation of Hi , H 2 

H i -> Hl(Hi,H 2 ) » = 1,2 (2.25) 
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with unit Jacobian 

det (||) = 1 (2.26) 

gives also 

v ND = VH[ x VH' 2 (2.27) 

Although the geometrical power of the surfaces H\ , H2 in the Nambu mechanics frame- 
work {D = 0) is lost, still there are interesting cases to consider. First the generalized 
Hamiltonians are not conserved but we have 

Hi = VD- Vffi 

H 2 = VD-VH 2 (2.28) 

and for D we obtain the equations: 

D = (VD) 2 + {D,H 1) H 2 } (2.29) 

if the relative orientation of the surfaces is fixed of positive sign, we get {D, Hi, H 2 } > and 
D increasing for all times. In general this is not true. 

We also see that if D is a surface orthogonal to either Hi or H 2 then the latter are 
conserved. If D is orthogonal to both surfaces, then it must be parallel to V-f/i x X7H 2 and 
the orbit of the dynamical system (2.24) is again the intersection of Hi and H 2 surfaces which 
are both conserved. 

The decomposition of a general phase space flow vector field in its (ND) and (D) parts 
may help to understand qualitatively the phase space portrait behavior of the full system. 
Indeed, in principle, it is possible to solve explicitly for the phase-space trajectories of the 
ND part 

'x = v ND = V Hi x V H 2 (2.30) 
using the two integrals of motion 

H{x) = H t {x ) i = l,2 (2.31) 

By parametrizing their intersection we can specify the explicit time dependence of the motion. 

In ref.[l5] we reduced the evolution equation of the form (2.32) Nambu Mechanics in 
Hamiltonian-Poisson form as follows: 

x { = {x\Hi,} H2 (2.32) 

where the induced Poisson bracket 

{f,g}H 2 = e i ^ k d i f&gd k H 2 (2.33) 
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satisfies all the required properties like linearity, antisymmetry and the Jacobi identity: 

{{f,g}H 2 ,h} H2 +{{g,h} H2 ,f}H 2 + {{h,f} H2 ,g} H2 = (2.34) 

The Jacobi identity follows from the Fundamental Identity rel.(2.7). 

Since H 2 is conserved it defines a fixed surface £2 : H 2 = H 2 {x ) embedded in B? which 
can be considered as a two dimensional phase space. In certain cases we may find explicit 
local or global parametrization of H 2 with two parameters, (u, v) € A C R 2 : 

x i = x\u,v) (u,v)eA i = l,2,3 (2.35) 

such that (x l ) € £2 

H 2 {x\u,v)) = H 2 {x ) (2.36) 

The orbits on the £2 phase-space are determined by the Hamiltonian H\ : 

x l = {x\R x ^ 2 i = 1,2,3 (2.37) 

Comparing the relations on T, 2 



and 



we get: 



x = x u ii + x v v (2.38) 

'x = VFi x VH 2 (2.39) 

it ■ (x u x x v ) = VH 2 ■ d v H 1 (2.40) 

v ■ (x u x x v ) = -VH 2 ■ d u Hi (2.41) 

Since x u x x v and VH 2 are normal and collinear at every point of T, 2 using parameters u, v 
of the same orientation (outwards) we must have: 

x u x x v = p(u,v) VH 2 (2.42) 

for a positive function p(u,v) on £2, we thus get from rel. (2.40-2.42) 

u = — ^— d v Hi (2.43) 
p(u,v) 

v = d u H x (2.44) 
p{u,v) 

Relations (2.43 — 2.44) can be written as Hamilton-Poisson eqs on T, 2 : 

u = {u,#i} S2 (2.45) 

v = {^,#i}e 2 (2.46) 
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with the Poisson bracket on T, 2 : 

{/,<?}s 2 = —Uf (d u fd v g - d u gd v f) (2.47) 
p(u,v) 

We shall use instead the level-set Morse functions H 2 for the phase space surfaces, because 
it is not always easy to find appropriate surface parametrizations. 

We observe that from eq.(2.33), the Poisson structure of the H 2 phase-space, is described 
in the algebra of phase space coordinates x % : 

{x\x j } H2 = e^ k d k H 2 (2.48) 

which is preserved by the eqns. of motion of the ND part. As an example of Nambu- 
Hamiltonian dynamics and the reduced Hamilton-Poisson form, we will consider the Euler 
free Top US]. 

For the Euler free top there are two conserved quantities. The total angular momentum 
squared: 

H2 = \ (l\ + l\ + it) (2-49) 

and the energy: 

H ^2{i + T 2 + i) < 2 - 50 > 

where (h),i = 1,2,3 are the components of the angular momentum, in the body frame, and 
I±,i = 1,2,3 are the eigenvalues of the moments of inertia tensor. The Nambu-Hamilton 
eqns : 

U = {luH^Hz} = ii k d 3 H x d k H 2 (2.51) 
provide the correct eqs. of motion for the rigid body : 

h = (l-i)/^ (2.52) 

The 2d-phase space is determined by the initial conditions: 

ll(t) + ll{t) + ll{t) = 2H 2 = / 2 = /?, + / 2 2 + Z|, (2.53) 

which is an S 2 sphere in i? 3 . 

The induced Poisson structure of the coordinates in ii 3 is: 

{k,lj}s? = £ijkdkH 2 = e ijk l k i,j,k = 1,2,3 (2.54) 

the Poisson SO(3) algebras and eqs. (2.52) can be written 

ii = {l i ,H 1 } s2 . = 1,2,3 (2.55) 

In the following sections we are going to present a detailed investigation of the Lorenz 
and Rossler attr actors from the point of view of Dissipative Nambu-Hamiltonian Dynamics. 
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3 The Lorenz System and its Euler top Deforma- 
tion 



3.1 Lorenz Attractor from Dissipative Nambu Dynamics . 

The Lorenz model was invented as a three Fourier mode truncation (Gallerkin approximation 
(3] of the basic eqs. for heat convection in fluids in Reyleigh-Benard type of experiments 
|23[ |24"1 [25] The time evolution eqns. in the space of three Fourier modes x, y, z which we 
identify as phase-space are : 

x = a(y — x) 
y = x{r — z) — y 

z = xy — bz (3.1) 

where a is the Prandtl number, r is the relative Reynolds number and b the geometric aspect 
ratio . 

The standard values for a, b are a = 10, b = | with r taking values in 1 < r < oo. There 
are dramatic changes of the system as r passes through various critical values which follow 
the change of stability character of the three critical points of the system P\ : x = y = z = 
0, P± : x = y = ±y/b(r — 1), z = r — 1 

Lorenz discovered the non-periodic deterministic chaotic orbit for the value r = 28, which 
is today identified as a Strange Attractor with a Hausdorff dimension of {d = 2.06) [26]. 
Standard reference for an exhaustive numerical investigation of the Lorenz system is the 
book by Sparrow [27 , although there is an extended list of numerical and analytic work on 
the system. For a qualitative discussion of the physics we refer to [7] [Til IE]- 

We can have a first glimpse of the geometry of motion, for a volume element in flow under 
eq.(3.1), if we cast the Lorenz system in Matrix form as follows: 



d_ 
dt 



( x\ 



X 

y 

\ z J 



t 



\ 



\ 



= A 



-a a 
r —1 —x 
x -b 

and decompose the matrix A into its symmetric and antisymmetric parts A = A s + A 



J 



x 

y 

\ z J 



(x\ 

y 

V z ) 



(3.2) 



/ 



A s 



A n 



-a 

a+r 
2 





r—a 
2 



q+r 
2 



(3.3) 



a—T 
2 



X 



(3.4) 
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The matrix A s executes an expansion along an axis x', a contraction along the orthogonal 
y' axis, as well as a contraction along the z-axis. The antisymmetric matrix A as executes 
a rotation around z-axis with a constant angular velocity and a rotation around x-axis 
with an x dependent angular velocity 

The axes x',y' are rotated fixed axes which we find through a diagonalization of A s . We 
have thus the basic qualitative mechanism for chaos through stretching, twisting and folding 

mum ■ 

We turn our attention now to an important characteristic property of the Lorenz attractor 
system, i.e. its "Orbit localizability" within an attracting(trapping) ellipsoid. |7 
Consider the ellipsoid in E?\ 

S 1 = rx 2 + ay 2 + a(z - 2r) 2 (3.5) 

Along the orbits we have 

^ = f-VSi = -2a [rx 2 + y 2 + b(z - rf - br 2 } (3.6) 

If we consider fixed values of S\ such that the corresponding ellipsoid S\ = constant lies 
entirely outside the fixed ellipsoid 

S 2 : rx 2 + y 2 + b(z - r) 2 = br 2 (3.7) 

then we get 

f-VSi < (3.8) 

for all orbits intersecting Si, i.e. all of them are ingoing to the surface. We observe that S 2 
acts like an attractor region. There have been various attempts made to localize the Lorenz 
attractor, by convex surfaces, in order to get information about Hausdorff dimensions [26 j 
and other characteristics |24| [28] . 

The case of Rossler attractor has been studied also from this point of view [30 . 

We will deal with the issue of Localization and the existence of an attracting ellipsoid in 
our Matrix Model generalization of the Lorenz systems, in order to get attracting ellipsoids 
in higher dimensional phase spaces. 

Concerning analytical methods we observe that the system can be explicitly integrated 
out, through the use of Liouville operator. For any observable p(x,y,z) on i? 3 we get the 
time evolution on trajectories 

p = C-p (3.9) 

with 

C = a{y-x)d x + [x(r - z) - y]d x + (xy - bz)d z (3.10) 
The formal solution for the x-coordinates is (C Q = C(x Q , y ,z )) : 

x(t) = e tc °x (3.11) 
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Similarly for y and z respectively. Since the trajectories are smooth and bounded, we must 
have convergent Taylor series in time|25} EU 132] : 

oo ,k 

<t) = E^m ( 3 - 12 ) 

k=o K - 

with 

x k = |t=o fc = 0,l,2,... (3.13) 
at 

From rel(3.11) we get 

x k = C k Q x k = 0,1,2,- •■ (3.14) 

So we can get as many coefficients xt we like by simply looking at the linear in the derivatives 
d Xo ,d Vo ,d Zo , part of £*. 

We now proceed to describe the Lorenz system in the framework of section 2. The flow 
vector field v is analyzed into its dissipative and non-dissipative parts as follows: 

vd = (—&%,—y,—bz) = VD (3.15) 

with the "Dissipation" function 

D = -- {ax 2 + y 2 + bz 2 ) (3.16) 

and 

vnd = (cry,x(r - z),xy) = (0, y, z - r) X (-a?, 0, a) (3.17) 
From rel (3.17) we get the two Hamiltonians or Clebsch-Monge potentials H\,H2 

VHx x VH 2 = vnd (3.18) 

or 

1 r .2 i I „ „\2 



Hi = 2 [y 2 + {z-rY] (3.19) 



and 



x 2 



H 2 = az - — (3.20) 
The Lorenz system (3.1) can thus be written in the form ( r = (x, y, z)): 

'r = V-ffi x VH 2 + VD (3.21) 

In the Non-Dissipative part(ND) of the dynamical system 

f = VHi x VH 2 (3.22) 

the Hamiltonians H\, H 2 are conserved and their intersection defines the ND orbit. Moreover 
if we get the reduced Poisson structure (sect. 2) from H 2 we obtain the 2-dim phase space T, 2 
to be the family of parabolic cylinders with symmetry axis the y-axis: 

H 2 = constant = H 2 {f* ) (3.23) 
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S2 is thus given by 

z = z + (3.24) 
la 

with x , z D the initial condition for x, z. The induced Poisson algebra (rel. 2.33) is given by 

{x, y} H2 = d z H 2 = a 
{y, z} H2 = d x H 2 = -x 

{z,x} H2 = (3.25) 

The dynamics on the 2d-phase space £2 is given by Hi 

x = {x,Hi}h 2 
V = {v,Hi}h 2 

z = {z,H x } H2 (3.26) 

and Hi is an anharmonic oscillator Hamiltonian with (x/a,y) conjugate canonical variables. 
Using rel. (3. 19-3.20) we get on £2: 

Hi = ^ [V 2 + \ (- 2 " « 2 ) 2 ] (3-27) 

with 

a 2 = x 2 Q - 2a(z Q - r) = -2H 2 + 2ar (3.28) 

where -j plays the role of the mass. Depending on the initial conditions we may have a 
sine; le well (a 2 < 0, H 2 > or) or a double well potential ( a 2 > 0, H 2 < ar ) respectively. The 
trajectories , the intersections of the two cylinders , H\ and H 2 with orthogonal symmetry 
axes (x,y) may either have one lobe left/right or may be running from the right to the 
left lobe. This is reminiscent of the topology structure of the orbits of the Lorenz chaotic 
attr actor. 




50 

Figure 1: Intersecting surfaces for Lorenz Attractor 
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In fig. 1 we present the Lorenz Intersecting surfaces for initial values xq = zq = 0,yo = 1- 
So the orbits of the non-dissipative system are the seeds around which and through the action 
of dissipation, the actual attractor is organized. In fig. 2 we present such an example for orbits 
of the nondissipative sector of the lorenz system superimposed with the full Lorenz attractor 
for the same initial conditions and standard values of the parameters (<r = 10,r = 28,i> = 



We may observe directly the action of dissipation. Similar observations have been made 
by Haken, in a different physics context and Nevir-Blender [33, 34 without use of the induced 
Poisson structure in B? rel.(3.25) and the ensuing interpretation [15] (rel.3.26) Sparrow has 
analyzed the integrable non-dissipative system, adding perturbative dissipative terms [27 . 




The full Lorenz system does not conserve H\ , H2 and there is a random motion of the 
two surfaces against each other. Their intersection is time varying. In effect at every moment 
the system jumps from periodic to periodic orbit of the non-dissipative sector. Moreover the 
motion of the non-dissipative system around the two lobes, either left or right, can now jump 
from time to time from one lobe to the other. 

We would like to mention at this point , another interesting feature of the Lorenz sys- 
tem. Due to the linear character of the dissipative part of the Lorenz attractor ,we are able 
to transform the full system to a volume preserving flow with time dependent generalized 
Hamiltonians ( or equivalently intersecting surfaces). 

This is done as follows. For a general Nambu-Hamiltonian system with linear dissipation: 



8/3). 




Figure 2: Nondissipative Orbit and Lorenz Attractor 




x 



y 



(V-ffi x ViT 2 )i - ax 
(VfTi x ViT 2 ) 2 - py 
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z = VHi x VH 2 ) 3 - iz 



(3.29) 



we define new variables (comoving frame) : 

x = e' at u , y = e~ /3t v , z = e'^w (3.30) 

Using the gradient operators with respect to q = (u, v, w) we can check that the eqs. of 
motion (3.29) become 

e-& + ^^{u,v,w) = V f ffixV f fl 2 (3.31) 

where 

Hi = H i {e- at u,e~ pt v,e" lt w) i = 1,2 (3.32) 

Eq.(3.31) implies that the (u,v,w) phase space volume elements are conserved. In effect, we 
have a new time variable r = . \ , e^ a+/3+7 ^ 



dt 



(u, v, w) = Vfffi x VjH 2 (3.33) 



The Lorenz attractor can be brought to the form of rel.(3.33). In terms of our Hamiltonian 
reduction it is a Hamiltonian system with a time dependent Hamiltonian Hi(e~ at u, e~ l v, e~ bt w) 
on a dynamical two dimensional phase space which is defined by the surface H2(e~ at u, e~ f v, e~ bt w) 
constant which moves in R 3 . In the comoving coordinates u, v, w the evolution of the system 
is volume preserving. A direct way to see the time dependence of Hi, H 2 and D, is through 
their time evolution eqs. 2.28-2.29. For the Lorenz attractor system we have 

Hi = \ {y 2 + {z-rf) 
H 2 = az - - 

D = - X - {ax 2 + y 2 + bz 2 ) (3.34) 
Their normal vectors are respectively: 



Vffi = (0,y,z-r) 
VH 2 = (-x,0,a) 

VD = (-ax,-y,-bz) (3.35) 



Finally we get: 



Hi = -y 2 - bz (z 
Ho = ax 2 — abz 



D = a 2 x 2 + y 2 + b 2 z 2 - (a + r) xy + (1 - b) xyz (3.36) 
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Since the Lorenz attractor is bounded in phase-space, the coordinates x,y,z have maxima and 
minima, which implies that the surfaces Hi,H2,S have maximum and minimum positions. 
So the Lorenz attractor is bounded by the volume enclosed by the maximum and minimum 
positions of these three surfaces. Estimates of the maximum and minimum values of the 
functions Hi , Hi have been determined using the constraint of the bounding Lorenz ellipsoid 
as well as other surfaces. 

Closing we would like to point out that our analysis of the Lorenz system applies equally 
to similar dissipative systems in i? 3 such as the Chen [29 and the Leipnik-Newton [T3] 
strange attractor constructions. They all share the same structure with the Lorenz system. 
It involves a nonlinear, nondissipative, integrable part which admits an intersecting surface 
Nambu parametrization superimposed to a linear dissipative sector. Last but not least it 
should be remarked that the R 3 geometrical picture presented above is directly generalizable 
to dissipative systems with R 2n+1 phase space dynamics. In complete analogy it is expected 
that 2n intersecting surfaces (Nambu Hamiltonians) generate the integrable sector of the 
system. 

3.2 A Strange Rigid Body Attractor as Dissipative Dynamics 
of Intersecting Paraboloid with a Cylinder 

We notice, in passing, that we can obtain a dissipative Euler top deformation of the Lorenz 
system, through the addition of terms rel.(2.55) : 

r = {l\H x } S 2 + d l D z = 1,2,3 (3.37) 

Indeed as we can see, the Lorenz attractor system 

l\r-l 3 )-l 2 

lH 2 - bl 3 (3.38) 

can be written as a dissipative Euler Top through the choice 1% = I3 = I for any I > and 
ii = jj^j but only for the special relation a = r of its free parameters a, r. The corresponding 
Dissipation function D is given by: 

D = -al l l 2 + X - [ail 1 ) 2 + (I 2 ) 2 + (Z 3 ) 2 ] (3.39) 

The last condition a = r, comes about from the constraint that we add to the Euler top 
non-dissipative flow, vnd = (0, — Z X Z 3 , Z 1 ^ 2 ), an irrotational (gradient) flow : 

vd = (<t(Z 2 — Z^), tI 1 — I 2 , — 6/ 3 ) (3.40) 



1 1 = 

1 2 = 

1 3 = 
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with the condition, 

V x v D = (0,0,r-a) = (3.41) 

We performed some numerical experiments for the values of a = r. We find that there is a 
chaotic attractor for values a = r = (25,26,27) and b = 8/3. 

We introduce now, to a 2-parameter Euler top deformation of the Lorenz system 

h = (j ~ hh + °{h - h) 

h = (| - |) + rh - h 

h = (j - y) hh - bh ( 3 - 42 ) 

where we have chosen I3 > I2 > I\. For more general linear terms we obtain the Leipnik- 
Newton linear control system on the motion of a Rigid Body where for appropriate values of 
the coefficients a double Lorenz attractor is obtained [13 . 

Performing a similar decomposition in a conservative(non-dissipative) and dissipative 
sectors of our system we get 

vnd = ((^-^)hh + crh,(^-T) l ^ + r h,(T-T) 1 ^) ( 3 - 43 ) 

and 

v D = (-vh,-h,-bh) (3.44) 
we can find two Hamiltonians H\,H2 and the Dissipation function D: 

V.ffi x VH 2 = vnd 
VD = -(ah,l 2 ,bk) (3.45) 

In order to get the two characteristic lobes of the attractor already at this level we notice 
that we can write the non-dissipative flow of rel.(3.43) as : 

v ND = (0,/3l 2 ,jh-r)x(-h,Sl 2 ,a') (3.46) 

where 

* = Trh^ = Trh> s= T=w x (3 - 47) 



and 



- ^ (3.48) 



From rel(3.46) we get, 
an elliptic cylinder along the l\ axis and 



Hi = \ Wl + hk-r) 2 ] (3.49) 



H 2 = a'h -1 + 4 (3-50) 
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a paraboloid with elliptic cross section. The Dissipation function D is identical to the case of 
the Lorenz attractor: 

D = -\{all + l 2 2 + Ul) (3.51) 

From the chosen ordering 73 > I2 > 1%, we get (3 > 0, 5 < (= Lorenz case ) and a' > 0. 
The H2, an elliptic paraboloid is oriented towards the positive £3 axis. The intersection of the 
two surfaces H\,H2 accounts for the existence of the two lobes of the attractor. Depending 
on the initial conditions it is possible that these lobes are connected with trajectories of the 
non-dissipative system. In fig. 3 we present the intersecting surfaces for the Euler- Lorenz 
nondissipative system for the initial values xq = zq = 0, yo = 1 an d values of the parameters 
a = 10, r = 28, b = 8/3, h = 1, J 2 = 2, J 3 = 3. 




200 

Figure 3: Intersecting surfaces for Euler Top Attractor 

A detailed analysis with numerical experiments will be presented in a forthcoming pub- 
lication. In fig. 4 we present the nondissipative orbit for Euler-Lorenz system superimposed 
with its full strange attractor for the same set of initial conditions and parameters as in fig. 3. 
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Figure 4: Nondissipative orbit and Euler-Lorenz Attractor 

4 Rossler Attractor from Dissipative Dynamics of 
a Cylinder Intersecting with an Helicoid 

Rossler introduced a simpler than Lorenz's nonlinear ODE system with a 3d- phase space, 
in order to study in more detail the characteristics of chaos, which is motivated by simple 
chemical reactions [S]. 

The Rossler system is given by the evolution eqns: 

x = —y — z 
y = x + ay 

z = b + z(x-c) (4.1) 

with parameters a, b, c usually taking standard values a = b = 0.2, c = 5.1 or a = b = 0.1, c = 
14 for the appearance of the chaotic attractor. The two fixed points of the system are: 

_, , . , c ± Vc 2 — 4ab c± Vc 2 — 4ab c±Vc 2 — 4ab. . ... 

F±:(x±,y±,z±) = ( , , (4.2) 

2 la la 

and depending on the parameters a, b, c their stability character changes. For the standard 
values both are unstable saddle points with the F_ having an unstable 2d manifold whereas 
the F + an unstable Id manifold. Studying the Poincare maps in the x-plane and the bifur- 
cation diagrams for y n +i = fiVn) return map, one discovers the period doubling bifurcations 
to be the underlying mechanism for the road to chaos [1 Oj. 
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vnd = (~y - z - —,x,b) (4.4) 



z 

vd = { — ,ay,z(x-c)) (4.5) 



We turn now to the study of the Rossler system as a Dissipative Nambu-Hamiltonian 
dynamical system. 

The key difference with the Lorenz attractor is that the dynamics of the system is simpler. 
Chaos appears as random jumps outwards and inwards the single lob attractor. 

In order to get the three scalars, the two generalized Hamiltonians H±,H2 which are 
conserved and characterize the non-dissipative part and D the dissipation term : 

f = V-Hi x VH 2 + VD (4.3) 

we checked after some guess work, that we must subtract and add a new term in the first 
equation. Indeed we find for the two parts 

2 

t 
2 

where we impose the constraints 

V-v ND =0 (4.6) 

and 

Vxt? D = (4.7) 
So we must determine H\ , H2 and D such that 

ViZi x VH 2 = vnd (4.8) 

and 

VD = v D (4.9) 

For D we find easily 

D = \ [ay 2 + {x- c)z 2 } (4.10) 
To get H\,H2 we must integrate first the Non-dissipative system: 

t 
2 

y = x 

z = b (4.11) 

The general solution is : 

z 2 

x(t) = -6(1 + z(t)) + (x + 6(1 + z )) cos(t) - (y a + z Q + - b 2 ) sin (t) (4.12) 
z 2 (t) z 2 

y (t) = b 2 - z(t) + {y + z + -± - b 2 ) cos(t) + (x + 6 + bz ) sin (t) (4.13) 



z 

x = — y — z — 
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z(t) = bt + z Q 

To uncover Hi , H 2 we introduce the complex variable 



(4.14) 



with 



We obtain 



with 



w (t) = W!(t) +iw 2 (t) (4.15) 



wi(t) = x(t) + b(l + z(t)) 

w 2 (t) = y (t) + z(t) + Z -^-b 2 (4.16) 



w(t) = w -e u (4.17) 



Wo = w (t = 0) (4.18) 
We see that there are two constants of motion, the first one being : 

I w(t) \ = \w \ (4.19) 

and we define correspondingly, 

Hi = \ \w(t)\ 2 = l -{ x + h {l + z)f + l -{ y + z+ Z l-b 2 f (4.20) 
The second integral of motion is obtained through the phase 

w{t) =\w \ -e t(f ° -e lt =\w \ e^ (t) (4.21) 

or from (4.14) 

<p(t) - Z - = <p -^ (4.22) 
and we define appropriately the second constant surface H 2 

v + z + — — 6 2 

* = » ^ 1+K i 2 + ,) - z (4 - 23) 

We easily check that rel.(4.8) is satisfied: 

VH ± xVH 2 = (-y-z-^-,x,b) (4.24) 

The family of surfaces H\ and H 2 are a quadratic deformation of a cylinder and re- 
spectively a quadratic deformation of a right helicoid. Their intersection is the trajectory 
(4.12-4.14). 
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5 Quantization of Dissipative Nambu-Hamiltonian 
Dynamics 



The main motivation for us to introduce the framework of Nambu-Hamilton Dynamics into 
the Dissipative Dynamical systems, apart from the very useful general global topological 
information of the phase-space portrait, is the possibility provided by this framework, for the 
formal quantization of chaotic attr actors. 

Indeed, based on our recent work [T2],we are able to quantize systems, for which the two 
Hamiltonians H\,H 2 are general quadratic functions of the phase space coordinates. 

Hi = ^x i M ij x j 

H 2 = Ir'.X i,j = 1,2,3 (5.1) 

where M,N are real symmetric matrices. 

Since for the Flossier attractor case we determined the Nambu surfaces Hi , H 2 to be of the 
complex form rel. (4.20-4.23) of higher degree than quadratic (see 4.20-4.23 ), we will proceed 
to work out in detail, in the next and last section of this work, only the quantization of the 
Lorenz system . 

In the present section we shall work out the method of quantization for a general quadratic 
system. 

At the end of the section we will provide some practical way for a Matrix Quantization 
of general dissipative dynamical systems in i? 3 (see also ref. |17j). 

Our strategy will be the following. Once the decomposition of the total flow vector field 
between dissipative and non-dissipative parts has been established we have the following 
representations : 

dx -* -* -* 

— = V#i x VH 2 + VD (5.2) 

The dissipative part vd is physically determined by either an unspecified "external" envi- 
ronment or it corresponds to a phenomenological modeling of our ignorance of some internal 
degrees of freedom. 

We shall firstly quantize the non-dissipative system written in the form rel(3.26) 

x i = {x\H±} H2 (5.3) 

In [15] we proposed the quantization of the Nambu 3-bracket and the induced Poisson struc- 
ture in B? as follows: 

{x\xi} H2 = e^ k d k H 2 i,j,k = l,2,3 (5.4) 

Introducing a quantization(deformation) parameter % the phase-space coordinates x 1 go over 
to operators for X 1 , satisfying commutation relations of the form, 

[X\X j ] = ihe ljk P k (X) z, j = 1,2,3 (5.5) 
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For polynomial Poisson brackets in (5.4) the corresponding three operators P k , k = 1,2,3 are 
respectively polynomials. They are determined (non-uniquely) by the following requirements 
a ) The Jacobi Identity 

[X 1 ,P 1 ] + [X 2 ,P 2 ] + [X 3 ,P 3 ] = (5.6) 
(3 ) The Existence of a Classical Limit 

lim P k {x) = d k H 2 A; = 1,2, 3 (5.7) 

7 ) The Existence of a Casimir H2(h), {H 2 (Jfi) H2)-> 

[X\H 2 {K]\ =0 i = 1,2,3 (5.8) 

and the commutation relation (5.5) should permit the 

5 ) Diamond Property of Unique Ordering of Monomials 

(X*) m (Xi) n (X k ) 1 i,j,k = 1,2,3 m,n,l = 1,2,- ■■ (5.9) 

into a sum of terms with given fixed order: 

(I 1 )™ 1 (X 2 )™ 2 (X 3 ) ms , (5.10) 

This last property generalizes the Birkhoff-Witt theorem for polynomial enveloping algebras 
of Lie algebras. The algebras (5.5) are called Polynomial Lie Algebras |35| [39] . The case of 
quadratic H 2 leads after quantization to classical Lie algebras. So in case where: 

H 2 = \ X'N^Xi (5.11) 

we get 

[X\X j ] = ihe ijk N kl X l i,j,k,l = l,2,3 (5.12) 

as well as the Heisenberg type of quantum evolution equations in R 2 are 

ihx 1 = (5.13) 

where Hi is an ordered (say Weyl orderer) operator in X , X 2 , X s . Each ordering determines 
a unique quantization procedure. 

For example, for the quantization of the Euler top we get 

h 2 = \ + + (5.14) 
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and for the eqs. of motion: 

h = \{Y 2 -j){hh + hi2) 

k = l^-Y^ + hh) (5.15) 

with time preserved commutation relations [li, lj] = ifieijklk where li are the standard angular 
momentum QM operators. If we study the spin S system then the Us are (2s + 1) x (2s + 1) 
Hermitian matrices. 

We proceed, in what follows, to include the dissipative part in the quantization procedure 
in such a way, that we recover the correct classical limit when the quantization parameter h 
goes to zero. Quantum aspects of open dissipative systems have been the focus of research 
in many areas of physics [18] . A general operator quantization scheme for non-Hamiltonian 
systems was recently developed by Tarasov [T7j. In our present case for a quadratic "Dissi- 
pation function" D ( as it is the case for the Lorenz system) we have a unique way to recover 
the classical limit, by choosing 

D = ^SijXi (5.16) 

then 

ihX* = [X\Hi] + ihS ij X k (5.17) 
This form of quantization for quadratic Hi, as it happens with the Lorenz system 

Hi = '.V'.l/^.V' (5.18) 



gives the quantum evolution equation, 



X 1 = A) k l -{X i X k + X k X j ) + S ik X k (5.19) 

where 

A) k = i lm M j{ N km i,j,k,l,m = l,2,3 (5.20) 

with the correct classical limit. The Weyl symmetric ordering leads to a general recipe for 
the quantization of dynamical systems with polynomial flow vector fields, in phase space 
coordinates. We simply utilize operator phase space coordinates, symmetrize the right hand 
side of classical evolution eqs. and forget about the parameter h. Comparison of different 
orderings can be done using the commutation relations (5. 12). 

The algebra of coordinates in the case of the non-dissipative system is conserved in time 
under the unique time evolution of Hi, since the Casimir H2 commutes with any H\. Includ- 
ing dissipation, H2 and Hi fail to be conserved and the commutation relations rel(5.5) are 
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not preserved in time. In the next section we will consider this formal quantization of the 
Lorenz attractor and we shall present the emerging Matrix Lorenz Model. 

Concerning the physical meaning of this formal quantization program, we would like 
to stress the fact, that we are concerned here with quantum fluids, physical systems while 
macroscopic they keep their phase coherence(superfluids). The various quantum fluid motions 
are described by changes of action of the order of Ti. The same holds true for the energy 
dissipation mechanism at the same space-time scales. These conditions seem to hold in 
recently observed turbulence of superfiuid helium [ 4 He 3 He — B ]. Theoretical treatments 
of these experiments show that probably even the classical fluid turbulence originates at the 
quantum level [36], although discussion is going on about this point of view. We plan to 
come back on these interesting points in our future work. 

6 A Matrix Attractor for the Non- commutative 
Lorenz System 

The Quantization method for Dissipative Nambu-Hamiltonian mechanics proposed in the 
last section corresponds to the passage from Hamilton-Poisson mechanics to Heisenberg's 
quantum mechanical evolution eqs. or matrix mechanics [38 . Although the basic duality 
principle of Quantum mechanics particle <-> wave, is not directly transparent in Heisenberg's 
matrix mechanics the passage to the Schrodinger's picture, where the wave nature of quantum 
mechanical particles is automatic, gets easily established. In the case of Non-Dissipative 
Nambu-Hamiltonian system in R 3 we generalized Hamilton-Poisson mechanics by employing 
the induced Poisson structure on i? 3 by the 3-bracket : 

x = {x\H x ,H 2 } = {x\H x } H2 (6.1) 

The induced Poisson structure on B? is degenerate and reduces to a foliation of non-degenerate 
Poisson structures, by the values of the level set Morse function H 2 = constant. 

Our proposal is a Matrix equivalent quantization in the line of Heisenberg's Quantum 
Mechanics. The Lorenz system has both H\, H2 quadratic which is also the case for its Euler 
Top deformation as well as the Leipnik-Newton system (sec. 3.2). Indeed from section 3 for 
the ND Lorenz system 

A* = {x\Hi,H 2 } i = 1,2,3 (6.2) 

we have, 

H X = \ [y 2 + [z-rf] H 2 = *z-*. (6.3) 
The -f^-phase space is an infinite parabolic cylinder with symmetry axis along y-axis. Given 
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the initial conditions x ,y , z Q it is given by 



x 2 — x 2 



with minimum at x = 0, 



* + (6.4) 



x 2 

s < 6 ' 5 > 



The corresponding Poisson algebra of coordinates in R 3 is: 

{x,y} = a {y,z} = -x {z,x} = (6.6) 

Upon quantization we get for the X,Y,Z Hermitian operator coordinates of the quantum 
mechanical phase -space 

[X,Y] = iha [Y,Z] = -ifiX [Z,X] = (6.7) 

The X, Z operators can be simultaneously diagonalized and play the role of two position 
"coordinates" while Y is a conjugate "momentum" operator to X. 

Since H2 is, by construction, the Casimir of the algebra it must be real diagonal and 
for unitary irreducible representations, H2 must be a multiple of the identity operator. The 
evolution Hamiltonian has been chosen to be H± 

Hi = \ [Y 2 + (Z - r) 2 ] (6.8) 

The quantum evolution equations for the Non-Dissipative Lorenz system are 



ihx = 


[X,H 1 


ih 2 




zhY = 








ihz = 


[Z,H! 


}h 2 


(6.9) 



or using the algebra rel(6.7): 

X = oY 

: _ _XZ+ZX +f . 

~ XY+YX , 

Z = (6.10) 

We observe that H\ along with the Casimir H2 are conserved. 

This has the consequence that for any time t the operators Xt,Y, Zt satisfy the same 
commutation relations(6.7). 

The ND quantum Lorenz system is integrable as it is the classical one. It is a quantum 
quartic anharmonic oscillator. Indeed from rel.(6.4) we get 

(X 2 (t)-X 2 ) X 2 H 2 

X = oY = - G 2 XZ ^ ZX + raX (6.11) 
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The Hamiltonian H\ becomes 

Hi = ^ + ^{X 2 - 2{ar-H 2 )f (6.12) 
We can deduce the eqn. for X directly from first two eqns. rel. (6.10- 6.11) 

X = -^X 3 + (ar-H 2 )X (6.13) 

With initial conditions 

X(t = 0) = X X(t = 0) = oY (6.14) 

satisfying the algebra rel. (6. 7) 

Depending on the values of H2 we may have double or single well quartic potential 
(repulsive or attractive linear term) with a rich structure of physical states (discrete spectra, 
instantons,etc...). 

We go back now to the full Dissipative quantum system. According to our discussion in 
section 5 to have a natural classical limit we shall apply the recipe of adding the dissipative 
terms (which are linear) with the 1T1 prescription. 

The question of going into physical models of dissipation for the Lorenz system, in the 
framework of single mode lasers, has been discussed in [40, 4b. There, a specific form 
is assumed, of quantum dissipation which is phenomenological. The general question of 
dissipative quantum mechanics is discussed in several places [18]. Nevertheless there is no 
universally accepted way for improving Quantum Mechanics including the environment into 
the Schrodinger or Heisenberg picture(inclusion of the observer !). 

According to sec. 5 we get 

X = a(Y-X) 

A XZ+ZX ~ ~ 

Y = + rX -Y 

2 

Z = XY + YX _ b z (6.15) 

Since the Casimir H2 is not any more conserved, if we assume the commutation relations 
(6.7) for the initial conditions X ,Y ,Z they will not be preserved in time. Indeed the role 
of dissipation is exactly the opposite of QM which assumes the incompressibility of phase- 
space. 

We shall now make a motivated guess, that since the classical dissipative Lorenz system 
has contracting volumes the same occurs under certain conditions, to be determined, for the 
quantum case as well. Motivated by this guess, which says that the effective quantum phase 
space is compact ( like in the Quantum Euler Top case) we proceed to study the Operator 
valued(Quantum) Lorenz system of eqs. under the assumption that X,Y,Z are hermitian 
N x N matrices for appropriate N = 2,3,.... We make, in effect, the proposal to study 
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the Quantum mechanical system in the approximation of a finite dimensional Matrix model 
(Finite Quantum Mechanics [43, 44]). 

In section (3.2) we generalized the Lorenz system by adding some new terms, which makes 
the total system resembling as a deformation of the free Euler Top. 

Moreover we observed that the new system has similar attractor structure with the Lorenz 
one. The new attractor is bigger in size but quite similar to the Lorenz one. We propose to 
reverse the picture and consider the new system as a dissipative deformation of the integrable 
Euler Top system (see the Nambu-Hamilton representation in section 2) which breaks the 
rotational symmetry down to a discrete Z<i symmetry (x, y, z) — > (— x, —y, +z). 

From this point of view the quantization of these systems can be implemented, through the 
use of N x N Hermitian matrices where N = 2s + 1, where s is the spin of the Euler top ( see 
rel. 3.14-5.15 ). The quantum Lorenz system then is a special case for 1% = I3 = /, I\ = 
with added dissipative terms. 

It is easy to check that in this case the multidimensional ellipsoid 

Sf = tr(rX 2 + 0Y 2 + a(Z - 2r/) 2 ) (6.16) 

appropriately chosen to be completely outside the ellipsoid 

Sf = tr(rX 2 + Y 2 + b(Z - r/) 2 ) (6.17) 

is in fact an attractor with any orbit passing through getting trapped. 

The proposed Matrix Lorenz system for Hermitian N x N matrices X,Y,Z rel. (6. 16) 
gives a non-linear ODE system with 3iV 2 equations. We observed that there is a global 
U(N) symmetry group in the system. For any U G U(N), we have that X — » UXW, Y — ► 
XJYU\ Z -» UZrf will be a solution if X,Y,Z is one already. As a result the invariant 
observables in the Matrix ODE v-system are traces of monomials with the simplest to consider 
being of the form 

tr(X), tr(X 2 ), tr(Y), tr(Z), tr(XY) ■ ■ ■ (6.18) 

in terms of which the eigenvalues of X, Y, Z matrices can be expressed. Interesting plots 
appear to be the triples of eigenvalues (Aj, i = 1, • • • , N as functions of time : 

X(t) - [Xi(t)] i=1 ... >N 
Y(t) -» [//i(t)]i=i,...,iv 

Z(t) -» [vi(t)} i=1 ,.. >N (6.19) 

Preliminary numerical experiments for N = 2,3 show Lorenz attr actors with more interac- 
tions between the two lobes. On the other hand the Matrix Lorenz system (6.20) has an 
absorbing ellipsoid and so each volume element shrinks to zero asymptotically in time. The 
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interesting question is, what is the HausdorfT dimension of the Matrix attr actor and under 
what condition survives quantum mechanics or vice versa. 

Concerning questions of interpretation of the matrix model we point out that if X, Y, Z 
are diagonal at time t = they will remain so for any time. In this case we get N decoupled 
Lorenz systems. If we impose as initial conditions small off-diagonal entries, then we get 
weakly interacting Lorenz systems. On the contrary large values for off-diagonal terms are 
associated to strongly coupled Lorenz systems. There is also a hierarchy of patterns for initial 
conditions which is preserved by time evolution. We can have block-diagonal structures with 

N x N = (iVi x iVi) (JV 2 x N 2 ) © • • • (N k x N k ) (6.20) 

with small off-diagonal interactions in each block so there appears the possibility of having 
patterns of self-similar structures as N grows. The non-commutative phase space structure 
resembles the one of Matrix Theory [35]. 

Extensions of Lorenz systems has been also considered in the literature, with x,y complex 
numbers and z real. Various Lorenz-like non-linear systems were studied with weak dispersion 
and dissipation [13] 02] . 

The interaction between different Lorenz attractors with mismatch in their respective pa- 
rameters a, r, b has attracted the interest of engineering(electronics- telecommunication)biological 
and Chemistry scientific communities because of the phenomenon of spontaneous appearance 
of either synchronization or anti-synchronization in amplitude or in phase of the non-linear 
oscillator [35] . 

Networks with nodes Lorenz non-linear oscillators provide different universality classes of 
scale free networks [37]. We can extent in various ways our matrix model transforming for 
example the parameters cr,r,b into constant symmetric N x N matrices (diagonal or not). 

The question of introducing relativistic or not field theories by replacing the harmonic 
oscillators of free fields (and the Feynmann perturbative graphs) by chaotic non-linear oscil- 
lators has been proposed by C.Cvitanovic in connection to the chaotic field theory and with 
turbulence [35] . Our Lorenz Matrix Model is a concrete idea exploiting this direction. 

As far as the possible connection of the Lorenz Matrix model with the above growing 
interdisciplinary literature, we believe that it can serve as a new laboratory for advancing 
theoretical issues on chaos and turbulence both classical and quantum . 

7 Conclusions-Open Problems 

The main result of our present work is the demonstration of Dissipative Nambu-Hamiltonian 
mechanics as the conceptual framework that underlies specifically strange chaotic attractors 
both in their classical as well as quantum-noncommutative incarnation. Nambu's dissipative 
dynamics of intersecting surfaces reproduced the familiar well studied attractor dynamics 
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of Lorenz, Rossler and Leipnik-Newton in a very intuitive manner accounting of their gross 
topological aspects. Moreover we showed how it generates a novel deformation of the Lorenz 
attractor by embedding into it the Euler top dynamics. The detailed properties of our chaotic 
Quantum Euler Top Attractor will be investigated in the future. 

In all cases the dynamics of surfaces admits a noncommutative-quantum generalization 
through their fuzzification. The general idea of non-commutativity of configuration spaces 
for fluids, fields and of strings, provides a unifying framework for physics [19]. This, in effect, 
led us to argue for the existence and construction , in a systematic and consistent way, of 
a Matrix model for the quantum Lorenz attractor. It is a generalization of the well studied 
Lorenz attractor to higher dimensional phase space. 

This was essentially achieved through the decomposition of a general chaotic flow into its 
rotational (non-dissipative) and irrotational (dissipative) parts. It is conveniently described 
by three scalar functions : the so called Clebsch-Monge potentials (or Nambu Hamiltonians) 
H\, H2 as well as the Dissipation scalar function D. For each concrete physical system , which 
is described by such (3d ODE's) flows the identification of such a dissipative or "environment" 
term is provided. Generally this term, as was the case of the Rossler attractor, is not an purely 
irrotational vector field. It encompasses phase space volume preserving components which 
we have to subtract away and render it irrotational with the phase space portrait dynamics 
of the trajectories as transparent as possible. 

The non-dissipative term must exploit the basic topological characteristics of the attrac- 
tors via the geometrical intersection of the constant surfaces Hi,H2- 

Given the constraints above we could still have the freedom to determine the Dissipation 
surface D associated with the smallest possible energy loss. 

In such a splitting, the non-dissipative and integrable sector comprises the strongest term, 
an indicator of how close the entire system is to being integrable or equivalently how big is 
the violation of its integrability from pure dissipation. 

Unexpectedly it uncovers gross topological features of the full dynamical system (double 
lobe Butterfly for the Lorenz system and its Euler Top deformation) or single lobe for the 
Rossler attractor. 

The Quantum behavior of Strange attractors was built systematically through fuzzifying 
the classical intersecting surfaces of the ND sector . We demonstrated this for the simplest 
case of the Lorenz system with a linear dissipation together with the existence of an attracting 
fuzzy ellipsoid. We constructed a Matrix (hence noncommutative) Lorenz attractor as a 
dynamical system of many coexisting attractors in 3-dim. phase space, interacting with each 
other, through the off-diagonal terms of the coordinate hermitian N x N matrices. The 
physical picture that it implies is a concrete quantum mechanical model for the transition to 
chaos and turbulence. 

A decomposition into a real diagonal forms along with the diagonalization of unitary 
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matrices of the phase-space coordinate matrices will make explicit the role of the Quan- 
tum Mechanical phases and the wave nature of the system. Synchronization phenomena in 
the amplitude or phase for the coexisting attractors relate equivalently to important issues 
concerning decoherence phenomena due to dissipation. Many Questions arise concerning 
Hausdorff dimensions mutual information and their scaling dependence on N deserve to be 
dealt with in the future. Similarly the behavior of time evolution of the eigenvalue distri- 
butions of the coordinate matrices X,Y,Z are of importance for establishing the type of the 
generated chaotic ensembles. Obviously we still at the beginning of our investigation and we 
hope that this work just opens an interesting set of questions. 
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